Pair-wise decoherence in coupled spin qubit networks 
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Experiments involving phase coherent dynamics of networks of spins, such as echo experiments, 
will only work if decoherence can be suppressed. We show here, by analyzing the particular example 
of a crystalline network of Fes molecules, that most decoherence typically comes from pairwise 
interactions (particularly dipolar interactions) between the spins, which cause 'correlated errors'. 
However at very low T these are strongly suppressed. These results have important implications for 
the design of quantum information processing systems using electronic spins. 
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A worldwide effort is presently on to make nanoscale 
solid-state qubits, whose purity and reproducibility can 
easily be controlled. Microscopic spins, existing in molec- 
ular magnets quantum dots 0] and semiconduc- 
tors Q, or in doped fullerenes Q, are a leading can- 
didate for this. In some of these systems (notably mag- 
netic molecules) , the individual qubit properties are con- 
trolled by chemistry instead of by nano-engineering (the 
'bottom- up' approach with spin Hamiltonians and 
inter-molecular spin couplings known and controlled to 
at least 3 significant figures. Spin also possesses other 
advantages - information can be encoded in the topolog- 
ical spin phase, with no need to move electrons around. 
Using spins for quantum information will ultimately re- 
quire (i) detecting and manipulating single-spins, and (ii) 
understanding and controlling decoherence. 

Single spins have been detected in a few ingenious ex- 
periments [6j, but we don't yet have a general-purpose, 
single-spin detection/manipulation tool, analogous to 
single atom STM/AFM. Consider, however, an array of 
spins, each having a low-energy doublet of states whose 
splitting is easily controlled by a magnetic field. Even 
without addressing individual spins, one can still demon- 
strate coherent qubit operation, using external AC fields 
to promote resonant transitions between levels, and pulse 
sequences (e.g. spin echo) to manipulate the phase and 
measure decoherence rates. This approach is well known 
for room-temperature bulk NMR quantum computing 
Q. Here we treat the case of electronic spins which, 
unlike nuclei, can be highly polarized at low T. We in- 
troduce a formalism allowing the description of any set 
of spin qubits obtained by truncation to low energy of a 
larger system, showing how the low-T decoherence rate 
can be dramatically reduced, even for a network of mu- 
tually coupled qubits. To be specific, we treat the case 
where the qubit is obtained by taking an anisotropic high- 
spin nanomagnet [|| with easy axis z, subject to a large 
transverse field Hj^ (Fig. [J), giving a low-energy dou- 
blet of states with easily controllable energy separation, 
2A Q (H^) [Fig. [2(a)]. To make quantitative and testable 
predictions, we then calculate the spin-echo decay rate in 



a network of Fes molecules. This is a clean, crystalline 
and stoichiometric chemical compound [1, 0] , where the 
inter-qubit and the qubit-environment interactions are 
known accurately, and it should be a good 'benchmark' 
for quantitative test of the theory. 

(i) Effective Hamiltonian: In a transverse field the 
effective spin Hamiltonian of the Fes molecule, with total 
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FIG. 1: (Color online) In a strong transverse field Hx, the two 
potential wells of an easy axis spin system approach each other 
on the Bloch sphere, (a) The spin anisotropy energy for an Feg 
molecular spin with Hi along y, easy axis z and hard axis x, 
when jj, Hy = 2.5 T. The low-lying states \Z±), are approx- 
imately localized in the two potential wells. The quantum- 
mechanical eigenstates are symmetric and antisymmetric su- 
perpositions of \Z±) (see text), separated by the tunneling 
gap 2A . (b) The resonance experimental set-up and the spin 
states on the Bloch sphere; (c) At T <C A /&b only the lowest- 
energy eigenstate is populated, \S) = 2~ 1 ^ 2 (\Z+) + \Z-)). A 
short /xwave pulse prepares the system in the \Z+) state (tt/2 
rotation, corresponding to 1/4 of a Rabi oscillation). The 
spin then tunnels coherently between \Z+) and \Z~) at a fre- 
quency 2A /h. The effect of static inhomogeneities in A can 
be compensated by a 7r-pulse in a spin-echo sequence (not 
shown) . 
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FIG. 2: (Color online) Fes in a transverse field (a) Tun- 
neling splitting 2A of the lowest doublet, as a function of 



U t \ = n gl^S 2 / AirV c — 0.127 K between nearest neigh- 
bors (here V c is the volume of the unit cell). However 
we are interested in the effective interaction between 
the qubits - this acquires a field-dependent and highly 
anisotropic tensor form V lJ (Hx) = g'TJ^g? when writ- 
ten in the truncated qubit basis. There may also be ex- 
change interactions [l2| between the molecules - these 
however have never been observed in Fes. Thus our low- 
T Hamiltonian for the dipole-interacting molecules be- 
comes: 



iieimg bpiiLLing za D ui me luweaL uuuuiei, dja a iluicliuii ui i 

H y . (b) Effective g-factors for Fes as a function of the ratio H c s — — ^ ^ 2Aq(Hi)s^ H — ^ ^ ^ ^ 



A (Hy) /Q . In both figures, the values at fi H y = 2.5 T are 
indicated by a black arrow. Our decoherence rate calculations 
are valid when A ^> Ud, and do not apply in the grey areas 
shown. 



spin S = 10, is controlled by crystal and external fields: 
H(S) = -DSl + ES 2 X + Ki(Si + St) - 7 S • H ± , (1) 



The easy/hard axes are along z and x respectively; here 
D/k B = 0.23 K, E/k B = 0.094 K, K£ jk B = -3.28xl0~ 5 
K, 7 = g e [A B [A a nd q e « 2 is the isotropic g-factor of the 
spin-10 moment jlOj. Henceforth we orient the external 
field along y, ie., Hi — * H y ; this can tune the tunnel- 
ing splitting 2A Q between the two lowest eigenstates over 
seven orders of magnitude [Fig. Ufa)]. Higher excited 
states are separated from the lowest doublet by a large 
gap Q ss 5 K, and for k B T <C VL Q the giant spin can be 
truncated to an effective spin- 1/2. 

Given the two lowest doublet eigenstates, |<S(Hj_)) and 
|.4(Hj_)), we write the quasi- localized states in the min- 
ima of the energy potential as \Z±) = 2- 1 / 2 (|5) ± |.4)) 
[Fig. Ufa)]. Defining the states \X±) = 2- 1 / 2 (|S) ± i\A)), 
we see that \X±) and \Z±) have the maximum (positive or 
negative) spin expectation values along x or z within the 
2-dimensional subspace of the qubit, where |<S) and \A) 
are the basis states. All of these Fes spin states can be 
calculated by numerical diagonalization of (JTJ) . We then 
define define qubit spin-1/2 operators s x , s y , s z , such that 
s z \S) = 1/2\S), etc. 

To describe the magnetic moment of the truncated spin 
qubit we introduce an effective (/-tensor g, also operat- 
ing in the qubit subspace, and defined so that the qubit 
magnetic moment [ljj m*(H) = (J, B Ylv 9nv(&)s v . In 
the geometry studied here, with H along y, g is diago- 
nal, with components: 

g x = g e ((X + \S x \X+) - (X-\S W \X-)); (2a) 
~g y = g e ({S\S y \S) - (A\S y \A)): (2b) 
g z = g e {(Z+\S z \Z+) - (Z_\S Z \Z_)). (2c) 

Numerical evaluation of g shows it to be highly 
anisotropic and field-dependent [Fig. [2fb)]. 

Consider now two spin-10 Fes molecules. The standard 
dipolar interaction | J2i^j S M v U^S^Sl has strength 
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(3) 

where the spins form a triclinic lattice [9j , and we choose 
the qubit z axis of quantization to be along the field 
(ie., along the original y axis). The full Hamiltonian 
also includes local spin-phonon and hyperfine couplings, 
whose detailed form we will not need here. 

(ii) Decoherence: Most general discussions of decoher- 
ence in qubit systems concentrate on '1-qubit' decoher- 
ence, ie., that coming from the interactions of individual 
qubits with the environment 13, Id , IH . In an insulating 
magnetic system like Fes both nuclear spins and phonons 
will contribute to 1-qubit decoherence [l6j . However this 
is not the only possible decoherence source. In a multi- 
qubit system one can have 'correlated errors', from pair- 
wise qubit interactions. A few analyses of this have been 
done [171 ]: depending on what model is chosen, these in- 
dicate that when qubits couple to the same bath, corre- 
lated decoherence is very serious, and may prevent error 
correcting codes from operating. 

In the set-up imagined here, decoherence will show up 
in measurements of the dephasing time, T^. We define the 
dimensionless decoherence rate as 7 ^ = %jTiA. a (the 'co- 
herence Q-factor' is then ~ 7r/ 7? j). We start by con- 
sidering the role of dipolar interactions - these are impor- 
tant because they exist in all spin systems, and cause cor- 
related errors via pair-wise spin interactions, exciting in- 
ternal modes of the spin system (ie., no external environ- 
ment). There are two ways to look at their contribution 
to 70. First, as a dephasing from dipole-mediated pair- 
flip processes, which in a resonance or echo experiment 
gives a homogenous absorption linewidth (Sui^) = T% ■ 
Second, as a scattering of the uniform spin precession 
mode off thermal magnons. In what follows we will as- 
sume that f7d/2A <c 1, ie., the dipolar interaction <C the 
tunneling splitting - only then will dipolar decoherence be 
small enough to make an experiment worthwhile (in Fes 
at, eg., HoHy = 2.5 T we have £/d/2A « 0.16). Since 
any reliable measurement of the decoherence rate will 
involve time-resolved spin-echo experiments, we neglect 
static inhomogeneous broadening in the calculations. 

If U d /2A <C 1 and k B T > A D , the contribution ^ 
due to pair-flip processes can be expressed in terms of 
the second moment of the homogenous absorption line, 



3 



{SwD « i 

Vleck analysis 



by incorporating the ^-factors in the van 



vV\2 



(t7) 



1 - tanli 



A 




yy (2g e S? 



(4a) 



(g 2 x ~g 2 z )TZU (4b) 



with W*,, 



Vc(|r^'| 2 (5 Mi , -3r^r^)/|r y '| 5 . The next term, 
neglected in Eq. (j4"a|) , is ~ 0(Ud/ A D ). This approach 
fails when k^T <C A , since the line assumes a Lorentzian 
shape, where is not related to ($co 2 )- 

Consider now magnon-mediated dipolar decoherence. 
In the experimental set-up we imagine here, a resonant 
tipping pulse applied to the spin ensemble causes a subse- 
quent uniform spin precession (ie., coherent tunneling of 
the spins in the crystal between states \Z+) and l-ZL)). 
This is equivalent, in our effective spin language, to a 
magnon with wave vector q = 0, with gapped energy 
huj ~ 2A D when A D 3> Ud- However the dipolar in- 
teraction couples this magnon to other magnons [l9| . 
The magnon spectrum w q is calculated using standard 
Holstein-Primakoff transformations 20] applied to Eq. 
©. The lowest-order processes conserving both energy 
and momentum here are 4-magnon processes [cf. Fig. 
Efb)]; for k-gT <C Ud these are the only ones that con- 
tribute significantly (6-th and higher-order processes are 
~ 0(k-Q,T /Ud) 2 relative to these). Selecting the Ti terms 
where the total spin polarization is unchanged, we de- 
rive a magnon contribution 7? to 74, given (again, for 
A > U d ) by: 
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r qq'| 2 -^[^q] <5 (cj D + CJ q - CJ q < -W q _ q ') (5) 



Here the relevant 4-magnon matrix element is T^, = 

(l/4A0[£( 4 )(q,q') + /C< 4 >(0, q')], where /C( 4 )(q,q') = 
2K yy (q - q') - K zz (q) - K xx (q) and 
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FIG. 3: Feynman self-energy graphs for the q = uniform 
precession mode, interacting with magnon excitations. We 
depict (a) a 3-magnon process, which is ruled out by energy 
conservation, (b) a 4-magnon process, and (c) a 6-magnon 
process. 



^[nj is the usual Bose statistical weighting of the 
magnon thermal occupation numbers n q , and fto; q = 
(A q -4|S q | 2 ) 1 /2 W ith A q = 2A o -\[2K yy {0)~K zz {q)- 
K xx (q)] and B q = \[K zz {q) - K xx (q) - 2iK xz {q)] for 
any q. The magnon analysis requires spin polarization 
close to unity, ie., k^T < 2A D , and therefore comple- 
ments the 'van Vleck' approach from the low-T side. 
It's interesting to notice that the two methods yield dif- 
ferent T-dependencies, ~ exp(— A /k^,T) for 7^ v and 
~ exp(— 2A /kBT) for 7™. The crossover occurs around 
ksT ~ A , where the lineshape turns from Gaussian 
(fceT > A ) to Lorentzian (k^T < A ). Full quantita- 
tive details will appear in a long paper |2~0 ] . 

To these 2-qubit decoherence processes one must also 
add 1-qubit decoherence processes to find the total 7^ 
that would be measured in an experiment. The contri- 
butions from interaction with phonons and nuclear spins 
have been calculated for Fes elsewhere (lijj . In large 
transverse fields nuclear spins give a rate 7^ s = i? 2 /2A 2 , 
where E a is the half-width of the Gaussian multiplet of 
nuclear-spin states coupled to the qubit. The phonon 
decoherence rate is given for Fes by [22j: 
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FIG. 4: (Color online) Dimensionless decoherence rates 7^ = 
h/T2A a as a function of tunneling gap 2A in Feg, at the in- 
dicated T. Thin lines: 7^ v arising from pair-flip processes, 
Eq. |gaj. We omit 7^ v at T = 0.05 K < U d /k B . Thick lines: 
7™ from magnon scattering, Eq. ©. The gap between the 
7^ v and 7™ lines is the crossover region between the valid- 
ity of the two methods. The dashed and dotted lines show 
respectively the phonon [7^ h , Eq. (7}] and nuclear (7^ S ) deco- 
herence rates at T = 0.05 K. The arrow indicates the optimal 
operation point of the Fes spin qubit at T — 0.05 K. Inset: 
7^ v and 7™ for an isotropic spin-1/2 on the same Fes lattice, 
as a function of the Zeeman gap g e fiBfJ- H y . 
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where M 2 AS (H y ) » ±D 2 \(A\S y S z + S z S y \S)\ 2 , with den- 
sity p — 1920 kg/m 3 and sound velocity c s = 1386 m/s. 

Fig. H] summarizes the results of all these calculations. 
Except at very low T and large A Q , dipolar decoherence 
completely dominates over nuclear and phonon decoher- 
ence. The optimal operating point at T = 0.05 K, having 
minimum total decoherence, is found where 2A Q « 1 K 
« 20 GHz, at fi H y rj 2.6 T; here the decoherence qual- 
ity factor is Q$ ~ 10 7 , with a coherence time T2 ~ 1 ms. 
Since 7^ h oc coth(A Q /fcBT) is essentially T-independent 
in this regime, while dipolar contributions to 7^ still vary 
strongly, a further decrease in temperature would not 
substantially decrease 7^, but would allow operation at 
lower frequencies. It is instructive to compare the results 
for the Fes qubit, with the case of a fictitious isotropic 
spin- 1/2 on the same lattice (inset Fig. SJ, obtained 
by setting g x = g y = g z = g e and replacing 2A D by 
g e fj,BfJ- H y in all our formulas. We then deal with simple 
Larmor precession instead of spin tunneling. There is a 
clear reduction in 7^, but no trivial proportionality fac- 
tor, because of the strong variation of g^(H y ) in the Fes 
qubit [Fig. Mh)}- 

We now consider the more general implications of these 
results. Note first that even for the set-up considered 
here, correlated errors cause large decoherence - we sus- 
pect they will even more strongly affect higher-order en- 
tanglement between the qubits. This is in line with the 
results in the quantum information literature [Tt| . but 
the concrete calculation here reveals a surprising feature, 
viz., that at very low T, the contribution of correlated 
errors can be made much smaller then the single-qubit 
errors coming from hyperfine and spin- phonon couplings. 

Dipolar interactions are dangerous for spin qubit de- 
sign, but they are hard to screen. One way to reduce their 
effects (apart from going to very low T [23j) would be to 
go to lower dimensional spin networks; recent progress in 
attaching and assembling nanomagnets on surfaces 24j | . 
or even in chain structures [25| might then yield viable 
architectures. Designs in which dipolar inter-qubit cou- 
plings can be made small - eg., low-spin systems like V15 
or Ci^Ni [2^| - and where the inter-qubit couplings re- 
sponsible for information manipulation can be switched 
on and off, are clearly favored. 

We believe that we have captured the intrinsic deco- 
herence processes in networks of coupled spin qubits, ex- 
tending to the case where the qubit is the low-energy 
truncation of a larger system. Spin-echo experiments on 
well-characterized systems like Fes would give a stringent 
test of the theory. Perhaps more important, such exper- 
iments would allow exploration of different spin network 
architectures, even before the manipulation of individual 
spins in such networks becomes possible. 

We thank A. Burin, W. N. Hardy, A. Hines, A. J. 
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